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Abstract 

We study mass-transport models with multiple-chipping processes. The rates of these processes are 
dependent on the chip size and mass of the fragmenting site. In this context, we consider fc-chip moves 
(where k = 1, 2, 3, ....); and combinations of 1-chip, 2-chip and 3-chip moves. The corresponding mean- 
field (MF) equations are solved to obtain the steady-state probability distributions, P(m) vs. m. We also 
undertake Monte Carlo (MC) simulations of these models. The MC results are in excellent agreement with 
the corresponding MF results, demonstrating that MF theory is exact for these models. 
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I. INTRODUCTION 



The formation of structures at the nanoscale has attracted a lot of attention QjJ, 12D - The morpho- 
logical and statistical properties of these systems are governed by microscopic processes like ad- 
sorption and desorption, fragmentation, diffusion and aggregation, etc. Different non-equilibrium 
steady-states can be accessed and phase transitions induced if the rates of the these processes are 
varied. Several experimental and theoretical studies have focused on adatom and cluster diffusion 



[3 



However, there are relatively few studies which consider the fragmentation of clusters, 



expected to commonly occur during hyperthermal ion beam depositions and sputter depositions 



II 1 0l - 4 1 3fl . The bombardment of clusters by energetic ions induces large island boundary fluctua 



tions which cause multiple-chipping events [ k3,LL4|]. These systems have received limited attention 



to date because of the unavailability of experimental probes to observe fragmentation. Further, the 
occurrence of these events is often on a pico-second time scale, making collection of data difficult. 

The processes of fragmentation and aggregation significantly affect the growth mechanism by 
altering the number of small and large clusters. The steady-state cluster-size distribution in these 



systems is generally characterized by an exponentially decaying tail 11 



15|] . In contrast, systems 



with adatom and cluster diffusion exhibit steady-state distributions which are power laws 



[ay 



In many physical systems, the fragmentation kernel depends upon the masses of the chip size 
and the fragmenting cluster. For example, recent experiments on Au clusters sputtered from em- 
bedded Au nanoparticles report distinct chipping kernels for small and large clusters [16]. In 
a related context, groups, herds, schools and flocks of animals also exhibit size-dependent frag- 



mentatio n fl 1 VL 11811 . Similar observations have been made in the context of polymerization [19], 
gelation [19] and complex networks [2(1. In all these systems, the steady-state distributions are 
either exponentials or power laws or their combinations. 

Fragmentation is thus a ubiquitous phenomenon, observed in a variety of physical systems. In 
the present paper, we study conserved mass models with mass-dependent chipping - the fragmen- 
tation rates depend on the chip size and mass of the fragmenting site. We obtain the steady-state 
mass distributions of these models in the mean-field (MF) limit, and compare them with results 
from Monte Carlo (MC) simulations. The purpose of our study is to clarify the generic features of 
steady-state distributions in the presence of multiple-chipping processes. 

This paper is organized as follows. In Section 2, we study models with mass-dependent 
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fragmentation and aggregation processes. We focus on models with fc-chip moves (where k = 
1, 2, 3, ...); and combinations of 1-chip, 2-chip and 3-chip moves. In Section 3, we present MC 
results for these models, and compare them with the corresponding MF solutions. We conclude 
this paper with a summary and discussion in Section 4. 



II. TRANSPORT MODELS WITH MASS-DEPENDENT CHIPPING 

Consider a mass-transport model on a discrete lattice where there is no adsorption or desorption. 
To begin with, masses mj(0) = 0, 1, 2, etc. are placed randomly at each site i with an overall 
mass density p. The evolution of the system is defined by the fragmentation kernel g m (n), i.e., 
the rate for a mass n m) to chip from a site with mass m. The n-chip then deposits on a 
randomly-chosen nearest neighbor. The mass of the neighbor adds up, while that of the departure 
site decreases, with the total mass of the system remaining conserved. Typically, the steady-state 
mass distribution [P(m)] of sites with mass m in these models shows either exponential or power- 
law decay with m. 

We study the above model within a MF approximation which keeps track of the distribution of 
masses, ignoring correlations in the occupancy of adjacent sites. Although the MF theory suffers 
from this deficiency, our MC simulations in Sec. 3 show that it gives an accurate description of the 
above model, even in the 1 -dimensional case. Let P(m, t) denote the probability that a site has 
mass m at time t. In the MF limit, P(m, t) evolves as follows: 

^ m OO 7712 

P(m,t) = -P(m,t) ^2 9m(mi) - P(m,t) ^ P(m 2 ,t) ^ ftmW 



dt 



mi=l m,2=l mi=l 



+ P(m + m 1 ,t)g m+mi (m 1 ) 

mi=l 

m oo 

+ 7", P(m-mi,t) 22 - p ( m 2,^)fi'm 2 ('^i), rn > 1, (1) 



mi=l m2=mi 



j oo rri2 oo 

P(0,t) = -P(0,t) P(m 2 ,t)^2g m2 ( mi )+ P{mi,t)g mi {m x ). (2) 

m,2=l mi=l mi=l 



dt 



Equations ([T])-© enumerate all possible ways in which a site with mass m may change its mass. 
The first term on the right-hand- side (RHS) of Eq. (Q]) is the "loss" due to chipping, while the 
second term represents the loss due to transfer of mass from a neighbor chipping. The third and 
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fourth terms are the "gain" terms which represent the ways in which a site with mass greater 
(lesser) than m can lose (gain) the excess (deficit) to yield mass m. The terms of Eq. © can be 
interpreted similarly. The above equations satisfy the sum rule 

d °° 00 

-^P(m,t) = 0, or ^P(m,t) = l, (3) 

m=0 m=0 

as required. 

Next we consider a specific form of the mass-dependent chipping kernel, g m (n) = D{m)n~ a , 
where a ^ is a parameter. This kernel allows for chips of all sizes, with small masses being 
more likely to fragment from a site than large masses. It is especially relevant in the context of 
fragmentation of sputtered clusters. As revealed by Auger and thermal-desorption spectroscopy 
measurements performed on these nanostructures, fragments of a few atoms display a large mo- 



bility on the surface which rapidly decreases with increasing cluster size lllll. 

Multiple-chip models may be interpreted as limiting cases of the above model. For a = 00, 
this kernel only gives rise to 1-chip processes. For large values of a, we expect the fragmentation 
to be dominated by 1-chip and 2-chip processes with different rates. As the value of a is reduced, 
the possibility of higher-chip processes becomes appreciable. Finally, in the limit a = 0, chips of 
any size are equally likely. Subsequently, we study 1-chip models, (1 + 2)-chip models, etc. as 
approximations to the model with g m {n) = D{m)n~ a . It should be noted that the n~ a -chipping 
kernel, unlike the multiple-chip models, includes the diffusion move, i.e., the movement of the 
entire mass m at a site. The steady-state distribution for this model will be discussed in Section 
III. 



A. fc-chip models and the case g m (n) = g(n) 

First, we study A;-chip models which act as building blocks to understand models in which chips 
of different sizes are allowed. The chipping kernel has the form 

g m (ri) = w5 n)k . (4) 

The corresponding rate equations for P(m, t), obtained by substituting Eq. © in Eqs. CQ)-©> are 
as follows (absorbing w into the time t): 

^-P(m, t) = -(1 + s k )P(m, t) + P(m + k,t) + s k P(m - k,t), m>k, (5) 
at 

—P(m,t) = -s k P(m,t) + P(m + k,t), m < k. (6) 
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Here, s k (t) = J2m=k P{ m i t) is the probability of sites having mass k or more. 

In order to obtain the steady-state solution P(m), the generating function Q(z,t) = 
Ylm=i zm P{ m i t) is computed from the above equations, and dQ/dt is set to in the steady 
state. P{m) is then obtained by evaluating the integral 

P( m ) = _L f dz 9!A m > i (7) 

where the contour C encircles the origin in the complex plane. The steady-state generating func- 
tion of the A;-chip model is 

n<~\ - z ( Sl ~ S ^ + Z ^ S2 ~ s s) + • • • • +z k s k (l - sQ 

Q{z) ~ (T^F) " (8) 

The conservation of mass requires that J2m=i ™>P{™) = P ? where p is the mass density. Putting 
dQ/dz\ = p, we obtain 

si + s 2 + + s k 

P = : • (9) 

1 - s k 



The corresponding P(m), evaluated from Eq. ©, comprises of k branches Ii22ll : 



fe-i 

p(m) = - ^+i)4 m_l)/ ^mod(m ) fc),r ( 10 ) 

8=0 

(We have defined s = 1 in the above equation.) Thus all the branches decay exponentially. The 
occupation probability of the i tn branch (i = ->■ k - 1) is Si = P(i) + P{i + k)+ P{i + 2k) + .... 
Notice that because of the nature of the A;-chip moves, the initial population Si in each of the 
branches is conserved at all times. Therefore, there are k conserved quantities in addition to the 
conserved mass. These enable us to determine the s/s in terms of p and the Si's. 

We briefly discuss the solution of the 1-chip model because of its ubiquitous nature. On substi- 



tuting k = 1 in Eq. (flOl) . the steady-state distribution simplifies to |23, 



24|]. 



P(m) =s?- s™ +1 = jL- (j^- p ) = aib?. (11) 

This 1-chip solution is actually valid for a wide range of mass transport models. It arises 
whenever the fragmentation kernel g m (n) = g(n), i.e., the chipping rate is independent of 



the mass of the departure site Q22Q . This can be verified by substituting Eq. (fTTT) in Eqs. (QQ)- 
©. Further, the above kernels can be written as a product of two non-negative functions, i.e., 
gm( n ) = f(n)h(m — n)/h(m) where h(x) is a constant here. They therefore satisfy the necessary 
and sufficient condition for the steady-state distributions to become factorizable. Evans et al. have 
shown that mean field theory is exact for this class of models |25|]. 



B. (l+2)-chip model 



We next consider a combination of 1-chip and 2-chip processes. Our MC simulations show that 
these mimic the g m (n) = D{m)n~ a model for a > 4. We consider the general (l+2)-chip model 
with the chipping kernel: 

g m (n) = wi8 njl 5 mjl + (w 2 S n>1 + w 3 5 nj2 )9(m - 2), (12) 

where wi, w 2 , w 3 are the respective chipping rates and 9(x) = 1 for x ^ 0. Notice that the above 
kernel has an explicit m-dependence, i.e., the 1-chip solution is not a steady-state solution for the 
corresponding rate equations except in special cases. 

Replacing g m (n) in Eqs. ©-©, we obtain the following rate equations: 

^-P(m,t) = -[wx{sx- s 2 ) + (w 2 + w 3 )(l + s 2 )]P{m,t)+w 2 P{m + l,t) 
at 

+w 3 P(m + 2,t) + - s 2 ) + w 2 s 2 } P{m - 1, t) 

+w 3 s 2 P(m-2,t), m>2, (13) 

j t P(l,t) = -[w 1 (l + s 1 -s 2 ) + (w 2 + w 3 )s 2 ]P(l,t)+w 2 P(2,t)+w 3 P(3,t) 

+ [w 1 (s 1 ~s 2 ) + w 2 s 2 }P(0,t), (14) 

|-P(0,t) = -[w 1 (s 1 -s 2 ) + (w 2 + w 3 )s 2 }P(0,t)+w 1 P(l,t) + w 3 P(2,t). (15) 
at 

As usual, we are interested in the steady-state solution of this model. The steady-state generat- 
ing function can be obtained using Eqs. (fl~3l- (fT5l) . and some algebra yields 

QM = (.6) 

N(z) = w 3 s 2 (l - si)^ 3 + [w 2 sx(l - s 2 ) - W!Sx(si - s 2 ) + w 3 sx(l - s 2 )]z 2 

+w 3 (sx-s 2 )z, (17) 
D(z) = -w 3 s 2 z 3 - [w 2 s 2 + wi(sx - s 2 ) + w 3 s 2 ] z 2 + (w 2 + w 3 )z + w 3 . (18) 

Further, the relation between the mass density p and sx, s 2 is obtained from dQ(z)/dz\ =1 = p as 
follows: 

w 2 s 1 + 2w 3 (s 1 + s 2 ) ^ 

P w 2 (l - s 2 ) + 2w 3 (l - s 2 ) - Wx(sx - s 2 ) ' 
To obtain P{m), we need to invert Q(z) using Eq. ©. This integral is done by finding the roots 

of the numerator N(z) (which is of the form zx quadratic in z) and the denominator D(z) (which 

is cubic in z), and rewriting Q(z) in the form of partial fractions. This procedure often proves to 



be very cumbersome. An easier route is to solve the recurrence relation for P(m), which can be 
obtained directly from Eqs. (fT3l- (fT5l) by setting the LHS to zero in the steady state. It may also be 
obtained by comparing the coefficients of the z m - terms on either side of Eq. (fT6l) . Some algebra 
yields the recurrence relation: 

p/ \ (w 2 + w 3 ) , s 2 (w 2 + w 3 ) +twi(si - s 2 ) . s 

Pirn) = P{m — 1H P(m — 2) 

w 3 w 3 

+s 2 P(m-3), m>3. (20) 
To obtain the solution, we assume that 

P{m) = Ax m . (21) 

Substituting this form in Eq. (|20l) results in a cubic equation for x: 

w 2 s 2 + w 1 (s 1 - s 2 



w 3 



s 2 + 



x-s 2 = 0. (22) 



w 3 

The roots of this cubic equation are denoted as x±, x 2 , x 3 , and are all real in this case Il26ll . (For 
brevity, we do not present their explicit forms.) 
The steady- state solution may then be written as 

P(m) = A x x™ + A 2 x™ + A 3 x™. (23) 



Notice that the recurrence relation in Eq. (1201) is valid for m > 3. The following choice of the 
coefficients A x , A 2 and A 3 ensures that P(0), P(l) and P{2) also follow Eq. d23l : 



(1 - s x )x 2 x 3 - (si - s 2 )(x 2 + x 3 ) + (s 2 - S3) 
(xi - X 2 ){Xi - x 3 ) 

A ^ _ -(1 - gl)xiX 3 + (Si - S 2 )(X! + X 3 ) - (s 2 - S3) ^ 

{xi - x 2 ){x 2 - x 3 ) 
A _ (1 - si)a;iX2 ~ (fi ~ s 2 )(x 1 + x 2 ) + (s 2 - s 3 ) 
3 (a?! - x 3 ){x 2 - x 3 ) 

Of the three roots, we find that | xi |, | x 2 |< 1 and | x 3 |> 1. In order to ensure that P(m) in 
Eq. (|23l) is meaningful, we set A 3 = 0. This condition, together with Eq. (fT9b . can be solved self- 
consistently to obtain si and s 2 in terms of w\, w 2 , w 3 and p. The coefficients A\, A 2 are thereby 
determined, and the steady-state distribution of the (1 + 2)-chip model is thus a combination of 
two power-law (exponential) functions: 

P{m) = A x x™ + A 2 x™. (27) 



As usual, the large-m behavior is dominated by the slower of the two exponentials. 

An alternative approach is to obtain A\ and A 2 in terms of {x\, x 2 , x 3 , S\, s 2 ) from the equations: 



P(0) = l-s 1 = A 1 + A 2 , (28) 
P(l) = si - s 2 = A x xi + A 2 x 2 . (29) 

We can then treat (x±, x 2 , x 3 , Si, s 2 ) as unknowns to be determined by 5 coupled equations as 
follows: 

(a) The Xi's must satisfy the cubic equation (|22l) . 

(b) The normalization condition provides the constraint 

1 — X\ 1 — x 2 

(c) Equation (fT9l for p. 

These equations can be solved numerically to obtain the unknowns in terms of wi, w 2 , w 3 , p and 
thereby the solution. 

It is useful to consider limits where the double-exponential function in Eq. (1271) reverts to a 
single-exponential form. 

(i) w 3 = 

If we substitute w 3 = in Eq. (fT2l) . the chipping kernel simplifies to 

flf m (n) = iMn,i<!>m,i + w 2 6 n> i8(m - 2). (31) 

This model has only 1-chip processes, but the 1-chip rate is different for sites with one unit of 
mass (m = 1) and those with two or more units of mass (m > 2). In this limit, the recurrence 
relation in Eq. (|20l ) reduces to 



P(m) 



s 2 H (si - s 2; 

w 2 



P(m-l), m > 2. (32) 



The resultant steady-state distribution has the following form: 



where 



a 2 b^ , m ^ 1, 

P(m) = < 2 ' (33) 

1 — Si, m = 0, 



W2{Si-S 2 ) _ Wx{Si - 8 2 ) + W 2 S 2 

a 2 — — — , o 2 — . (J4) 

wAsi - si) + w 2 s 2 w 2 
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We obtain si, s 2 as functions of wx,w 2 , p: 

_ -w 2 (l + p) + y/w%(l - p) 2 + Aw lW2 p 

2{Wx - W 2 ) 

S2 = ^ . (36) 

W 2 + [Wx - W 2 )Si 

The exponentially-decaying solution in Eq. (l33l) has a different slope from that in Eq. (fTT|) . In 
the limit iuj = u> 2 , the solution in Eq. (|33|) reduces to the 1-chip solution of Eq. (fTT|) . 
(ii) u>i = w 2 , w 3 arbitrary 

If we substitute W\ = w 2 in Eq. (fl~2l) . the chipping kernel reduces to 

g m (n) = (w 1 5 nil +w 3 5 ni2 )9(m-n). (37) 

This kernel has the functional form g m (n) = g(n), and the corresponding P(m) is the 1-chip 
solution in Eq. (fTTI) . 

C. (l+2+3)-chip model 

We now consider a model with the possibility of chipping 1, 2 or 3 units of mass: 

g m (n) = Wxdn^Srn^ + iw^x + w^n^dra^ + iw^x + w^n^ + Wedn^eim-^). (38) 

Substituting Eq. (|38l) in Eqs. ©-© yields the required rate equations, which we do not present 
here. Some algebra yields the generating function: 

Q(z) = §4' (39) 

N(z) = w 6 s 3 (l - sx)z 5 +[ - (w 2 + W 3 Sx)(s 2 - s 3 ) + w 4 (s 2 - s 3 ) 

+ (w 5 + w 6 )(s 2 - SxS^Z 4 + [-Wxsx{s\ - s 2 ) - (w 2 + w 3 )s 1 (s 2 - s 3 ) 

+ (W 2 - W 4 )(S 2 - S 3 ) + S1O4 + U7 5 + W 6 )]z 3 + [W 5 (sx - S 2 ) + Wq(s 2 - S3)] Z 2 

+w 6 (si - a 2 )z, (40) 
= -w 6 s 3 z 5 - [w 3 (s 2 - s 3 ) + (w 5 + w 6 )s3]^ 4 - [wx{sx - 8 2 ) + (w 2 + w 3 )(s 2 - S3) 
+ (wx + w 5 + w 6 )s 3 ]z 3 + (u>4 + w 5 + w e )z 2 + (w 5 + w e )z + Wq. (41) 
As before, the steady-state distribution is obtained from the recurrence relation for P(m): 
w 6 P(m) = —(w 5 + w e )P(m — 1) — (u>4 + w 5 + w 6 )P(m — 2) 

+ [wi(si - s 2 ) + (w 2 + w 3 )(s 2 - s 3 ) + (w 4 + w 5 + w 6 )s 3 ] P(m - 3) 

+ [w 3 (s 2 - s 3 ) + (w 5 + w 6 )s 3 ] P(m - 4) + s 3 w 6 P(m - 5), m ^ 5. (42) 



Substituting P{m) = Bx m in the above equation results in the following quintic equation: 



W 6 



+ I 3 

w 6 



— [si - s 2 ) H («2 - S3) + H S3 



H s 3 H (s 2 - s 3 J 

io 6 y w 6 



x — s 3 = 0. 



(43) 



This quintic equation cannot be solved explicitly in terms of radicals. However, we can generalize 
the alternative approach described for the (1 + 2)-chip model subsequent to Eq. (1271) . The resultant 
coupled equations (solved numerically) yield the required solution for given values of w\ to w e 
and p. Typically, P{m) for the (l+2+3)-chip model is a sum of three exponential functions. 

As before, it is instructive to examine some simple limits of this model, 
(i) w 3 ,w 5 , w 6 = 

This model has only 1-chip processes, but the chipping rates depend on the mass of the departure 
site: 



g m (n) = w 1 5 nA 5 mA + w 2 5 nA 5 rnt 2 + w 4 5 ntl 8(m-3). 



(44) 



The corresponding distribution can be obtained from the simplified version of Eq. (1421) . or the 
simplified Q(z) from Eq. (|39| ): 



P(m) = < 



a 3 bf, m > 2, 
si — s 2 , m = 1, 
1 — Si, m = 0, 



(45) 



where 



«3 



wl{s 2 - S3] 



wAsx - s 2 ) + w 2 (s 2 - s 3 ) + W4S3 



(46) 



[Wi(s! - S 2 ) + W 2 (S 2 - S3)] 2 + W4S3' " W4 

These results, along with the equation for p, enable us to determine s/s in terms of Wj's and 
p. It is straightforward to generalize Eqs. (f33])-(r34l) and Eqs. (H3T)- (|46l) to the 1-chip model with 

different rates for departure-site masses m — 1, m — 2, , m ^ k. 

(ii) w 6 = 

With this substitution, the chipping kernel of Eq. (|38l) reduces to a generalized version of the 
(1 + 2)-chip model discussed earlier: the 1-chip and 2-chip rates are now different for sites with 
m — 1, m — 2 and m ^ 3. The recurrence relation in this case also yields a cubic equation. The 
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steady-state probability distribution has the same form as Eq. d27l) . but the decay rates of the two 
exponential functions are distinct from those in Eq. (|27l) . In the limit w 2 = u>4 and w 3 = u> 5 , we 
recover Eq. (1271) . 
(iii) Wi = w 2 = u>4 and w 3 = w 5 

This case corresponds to the chipping rates being independent of the site mass. In this limit, the 
chipping kernel in Eq. (1381) becomes 

9m{n) = {wi5 n ^ 1 + w ? ,5 n: 2 + w^ n ^)d(m-n). (47) 

Hence, the steady-state distribution in this limit is the 1-chip solution of Eq. (fTTI) . 



III. MONTE CARLO SIMULATIONS 



In this section, we present Monte Carlo (MC) results for some of the models discussed earlier. 
All simulations were performed on l-d or 2-d lattices with periodic boundary conditions. The 
lattice sizes were L = 1024 in d = 1, and L 2 = 128 2 in d — 2. The initial condition for a 
run consists of a random distribution of masses with density p. We evolve the system using the 
chipping rate g m (n), and compute the mass distribution P{m). This quantity settles to equilibrium 
for t > 25000 MCS - we show results for P{m) vs. matt = 50000 MCS. The statistical data 
presented here was obtained as an average over 200 independent runs. 

First, we present results for the kernel g m (n) = g(n), discussed in Sec. 2.1. In Fig. 1, we plot 
P(m) vs. m obtained from \-d MC simulations with p = 5 and three different functional forms of 
9m(n) = 1,1/ 'n, n 3 e~°' 2n and p = 5. The MC data sets are numerically coincident with each other. 
They are also in excellent agreement with the 1-chip solution in Eq. ([141) (denoted as a solid line), 
which was obtained from the corresponding MF equations. Our subsequent results also show that 
the MC data is described very well by the solutions of the corresponding MF equations, even for 
d = 1. This demonstrates that the MF equations are exact in the present context 12511 . 

Second, we present results for the (l+2)-chip model discussed in Sec. 2.2. Fig. 2 shows 
the steady-state distributions obtained from l-d and 2-d MC simulations with p = 5, and W\ = 
6, w 2 = 0.5 and w 3 = 7. The solid line denotes the result in Eq. (ITTT ) - the corresponding roots 
and coefficients are specified in the caption. The oscillatory structure of P{m) arises as one of the 
roots is negative. For the wide range of parameter values we considered, one of the roots is always 
found to be negative. 
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Next, let us consider the arbitrary kernel g m (n) = D{m)/n a with D(m) determined from 
Y^n=i D{m)/n a = 1. Unlike our earlier models, this fragmentation kernel allows diffusion, which 
corresponds to the movement of the entire mass m at a site to a randomly chosen neighbor. In Fig. 
3, we show P{m) vs. m from MC simulations in d = 1,2 with a = 3 and p = 5. The solid line 
denotes the analytical result for the analogous (1 + 2 + 3)-chip model with Wi's specified in the 
caption. Clearly, the (1 + 2 + 3)-chip approximation captures the original model rather well - this 
is true for a > 2. 

Finally, in Fig. 4, we show P(m) vs. m from MC simulations in d = 1, 2 for the 1-chip model 
with different chip rates upto m ^ 5, i.e., g m (n) = 5 n ,i[wi5 TO) i + w 2 S m ,2 + ^3^m,3 + w A 5 m ^ + 
w 5 9(m — 5)]. We show results for p = 5 and the io/s are specified in the caption. The solid 
line denotes the analytical solution obtained by generalizing the solution in Eqs. (|45l)-(|46l). This 
solution is exponential for m ^ 4. 



IV. SUMMARY AND DISCUSSION 



We conclude this paper with a summary and discussion of the results presented here. We 
have studied the steady-state distributions [P{m) vs. m] for mass-transport models with multiple- 
chipping processes which depend upon the masses of the departure sites and the chips. These 
models are relevant for a variety of physical applications. In general, a site with mass m could 
have 1, 2, 3,.., m units of mass chip with different rates and aggregate with a neighbor. In this 
context, we study fc-chip processes (where k = 1,2,3,...) and combinations thereof. We undertake 
Monte Carlo (MC) simulations of these models, and analytically study the steady- state solutions 
of the corresponding mean-field (MF) equations. Our MC results are in excellent agreement with 
the MF solutions, demonstrating that MF theory is exact in this context. This is true even when 
the kernels are not factorizable Q25Q, and requires further investigation. 

The steady-state distribution of the A;-chip models has k branches, each of which decays expo- 
nentially with the same slope. We find that a large class of chipping kernels, where g m (n) is inde- 
pendent of m, gives rise to an exponentially-decaying distribution: P{m) = (l + p) _1 [p/(l + p)] m , 
where p is the mass density. This is also the MF solution for the 1-chip model, where one unit of 
mass fragments from a site and aggregates with a randomly-chosen nearest-neighbor. 

We have also studied models with a combination of 1-chip, 2-chip and 3-chip processes. The 
steady-state distributions of the (1 + 2)-chip model and (1 + 2 + 3)-chip model are sums of two 
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exponential functions and three exponential functions, respectively. We expect the steady-state 
distribution of the (1 + 2 + . . . + /c)-chip model to be a sum of k exponentially-decaying functions. 
We have also examined several limiting cases of the above models. Our conclusion is that the 
steady-state distribution is sensitive to slight changes in the chipping kernel. 

Finally, it is relevant to discuss physical processes which give rise to multi-exponential mass 
distributions vs. power-law distributions. Power laws have been observed in a class of mass- 
transport models where single-particle adsorption or chipping processes and (mass-independent) 
diffusion processes maintain a delicate balance between the lower and upper ends of the mass 
spectrum |7-^]. On the other hand, mass-dependent fragmentation and aggregation precludes this 
balance. Consequently, all our models exhibit exponential distributions. These would be relevant 
in the context of ion-beam and sputter-deposited nanostructures, animal group distributions, 
polymerization, gelation and complex networks. 
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FIG. 1: Steady-state probability distributions [P(m) vs. m] from d = 1 MC simulations with three different 
forms of g m (n). The data sets are plotted on a linear-logarithmic scale. The details of the MC simulations 
are provided in the text. The mass density is p = 5. The solid line denotes the 1-chip solution in Eq. (fTTT) . 
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FIG. 2: Plot of P(m) vs. m for the (l+2)-chip model, obtained from MC simulations in d = 1,2 with 
p = 5. We used g m {n) in Eq. (fT2l with u)i = 6, = 0.5 and W3 = 7. The solid line denotes the solution 
inEq. $2J$ with xi = 0.8427, x 2 = -0.6498, A\ = 0.1475, A 2 = 0.1021. 
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FIG. 3: Plot of P{m) vs. m, obtained from MC simulations in d = 1,2 with kernel g m (n) = D(m)/n a 
[where D{m)~ 1 = Y^=i n ~ a ] for a = 3 and p = 5. The solid line denotes the result for the analogous 
(l+2+3)-chip model with w x = 1, w 2 = 8/9, w 3 = 1/9, w 4 = 216/251, w 5 = 27/251 and u> 6 = 8/251. 
The solution is a sum of three exponentials with x\ = 0.8372, x 2 = —0.0594 — 0.1263i, X3 = —0.0594 + 
0.1263i, Ai = 0.1590, A 2 = 0.0118 + 0.0090i, A 3 = 0.0118 - 0.0090i 
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FIG. 4: Plot of P(m) vs. m for the 1-chip model with different chip rates upto m ^ 5. The symbols 
denote MC simulation results in d = 1,2 with p = 5 and w\ = 1, = 2, = 3, = 4, = 5. 
The solid line denotes the generalization of Eqs. (|43T)-(|46T). The solution decays exponentially for m ^ 4, 

P(m) = 0.5133 x (0.7421)" 1 . 
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